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ABSTRACT 

A rigorous mathematical proof is presented for a multiline 
representation of the equivalent width of a molecular band which 
consists in the general case of n overlapping spectral lines. The 
multiline representation includes a principal term and terms of 
minor significance. The principal term is the equivalent width of the 
molecular band consisting of the same n nonoverlapping spectral 
lines. The terms of minor significance take into consideration the 
overlapping of two, three and more spectral lines. They are small 
in case of the weak overlapping of spectral lines in the molecular 
band. The multiline representation can be easily generalized for an 
optically inhomogeneous gas media and holds true for combina- 
tions of molecular bands. If the band lines overlap weakly the 
standard formulation of line-by-line method becomes too labor- 
consuming. In this case the multiline representation permits line- 
by-line calculations to be performed more effectively. Other useful 
properties of the multiline representation are pointed out. 

1. INTRODUCTION 

As more complete and precise i nformation about parameters of 
spectral lines becomes available (Rothman et al., 1987), the line- 
by-line method is used as a reference method in computation of the 
total vibration - rotation molecular and absorptance. This method 
assumes that the band consists in the general case of n overlapping 
spectral lines each of which is characterized by intensity S p half- 
width Yj , frequency of line center v9 and line profile a(v^-v\i=\ y 
2,...,n. The equivalent bandwidth W(n), used for a description of 
total absorptance of the molecular band in optically inhomogeneous 
gas media is given by 


W(n) = jav 


1 - exp - 


n O 

XTjCX(Vj - V) 


( 1 ) 


where is optical depth at the central frequency of i-th line. 

Line-by-line method assumes numerical integrating of the 
oscillating function, the number of the oscillations in the general 
case being equal to the number of spectral lines composing the 
band. In this case the calculation time rises with increasing of the 
oscillation magnitude of the absorption coefficient at line center, as 
the line overlapping is reduced. Thus, in the case of weak overlap- 


ping of the lines in the band the line-by-line method becomes too 
labor-consuming and then a new approach is needed to calculate the 
total band absorptance. The following multiline representation can 
be taken as the basis of such an approach. 

2. THE MULTILINE REPRESENTATION 


The representation (Sakai et al., 1 964) is known for the equiva- 
lent width W(2) of the band consisting of two overlapping spectral 
lines: 

2 

W<1)£ W, + W 2 - J dvfl {l - ex P [-x i a( v° - v)]}, (2 ) 
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where 

W i = | dv { 1 “ ex p[“ T i a ( V P “ v )]}» 0) 


is the equivalent width of the i-th line. The relation (3) follows from 
the obvious identity 1 - exp( -qj - q 2 ) = =(l-exp(-qj))+(l-exp 
( — q 2 )) — ( 1 — cxp( — q j ))xx( 1 -exp(-q 2 )) w'hich takes place for any 
material number qj and q 2 . It may be proved by means of the 
mathematical induction method with respect to n, that for n > 2 we 
have 
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- X n( 1 “ ex p( _t li m )) + - + ( _ 1 ) ntl n(l-exp(-q i )). 
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Substituting in (1 ) and using q^^x^v^-v) we obtain the follow- 
ing representation for the equivalent bandwidth; 

W(n) = tw i Jdvfl{l-exp[-x lm a(vP m - v)]} — 
11 MJ2 m_1 (4) 
+ - + (-l ,+1 X Jd.nfl-expK xxa(vf-v)]} 

i] i 1 m 1 

+ ... + (-l) nl Javjqjl-exp^Tj xxoc(v° - v)]J. 
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The first term in the right-hand side of (4) represents the 
equivalent bandwidth W N (n) in the nonoverlapping line approxi- 
mation. In this approximation each line corresponds to its own 
ensemble of photons and these ensembles do not blend. The second 
term takes account of the overlapping of two lines or blending of the 
photons from two different ensembles, ets. By analogy with well- 
known problem of electrodynamic the proposed relation (4) can be 
referred to as multiline representation of the equivalent width of the 
band consisting in the general case of n overlapping spectral lines. 
In spite of its awkward form, multiline representation (4) has 
proved to be useful in practice for the following reasons. Represen- 
tation (4) is an asymptotic series expansion for W(n), from (1), in 
which the equivalent bandwidth W N (n) in the nonoverlapping line 
approximation is the principal term because the 1-th term in the 
right-hand side of (4) decreases as jdv(n t(v9 m-v) with increas- 
ing of sufficiently large distances between the line centers and 
decreasing of sufficiently small optical depths in the line centers 
(i.e. as the nonoverlapping line approximation becomes more 
valid). Accordingly, the estimating of W(n) by summarizing only 
few first terms of (4) will be more exact as the nonovcrlapping line 
approximation becomes more valid, i.e. as the line-by-line method 
becomes more labor- consuming. Besides, as long as the minor 
terms of (4) in this case will make insignificant contribution, as 
compared to the principal term W N (n), it is possible to estimate 
them based on physical meanings. 

3. APPROXIMATION FOR MINOR TERMS OF (4) 

Let us arrange the band lines in the order of decreasing optical 
depths Xj (X]>X2- use approximate relation (Ishov et al., 

1990): 

fdv J {l - exp[-TjOt(vP - vj]| = 

i -1 

(5) 
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Then it may be shown with the mathematical induction method 
with respect to n that instead of (4) we have: 


W(n) = W^(n) - A, A * 


X w i 


i l 





k i * 1 



(6) 


The approximate relation (5) is found to be useful in the case when 
radiative transfer occurs mainly at the frequencies of the spectral 
line with maximum x ] on a background of the comparatively 
smooth changing of the wings of other lines, included in the integral 
of (5), and for which X;<X,. Of course, one can construct other 
representations for A of (6). However, the above representation is 
preferable because in addition to elementary' functions it contains 


functions W- only, describing isolated line absorptance, the values 
of which have already been calculated with calculating the princi- 
pal term of multiline representation (4). Clearly, that other condi- 
tions being equal, relation (5) is worst satisfied when all X^j. 
Besides, with weak overlapping of spectral lines the overlapping 
occurs primarily in the line wings which according to the Voigt 
profile conception (i.e. given statistically independent Doppler and 
Collision mechanisms of the spectral line broading) have the 
Lorentzian character. Therefore we can limit oneself to the Lorentz 
profile case when testing the approximate relation (6). It should also 
be noted that for real vibration - rotation bands the variations of the 
half- widths ^ and the distances between neighbouring lines Iv9 +1 - 
v9l are negligible in comparison with those of the intensities 
Sj(i.e.X|). Therefore it seems reasonable to illustrate the results 
obtained using the model critical case rather than single examples 
of the specific bands. The case is a band, consisting of n identical 
equidistant Lorentzian lines, which at n-°o represents Elsasser 
model band. The testing results are shown in Fig. 1 . It illustrates that 
using the approximate formula (6) significantly increases the 
accuracy of W(n) calculations in comparison with W N (n). 
Farthcrmore, there is a “corridor” of values 6 and x, in which 
approximate representations (6) proves to be extremely accurate. 

4. PRINCIPAL TERM OF (4) 

ITie equivalent width of an isolated spectral line W(x) with 
optical depth x at the line center is expressed through dimensionless 
absorption function Mq(x) and its derivatives M k (x) with respect to 
X (k= 1 ,2,. . .) for I xirentz (L), Doppler (D) and Voigt (V) line profiles 
which are mathematically well-studied (Ivanov, 1973). High- 
precision calculation algorithms have been developed for 
and M k (x) functions (Hummer, 1981, Ishov, 1982). For 
functions we have power for X *0 and asymptotic for x~*oo (a+«>) 
series expansions (Ivanov, 1973;Hummer, 1982) and Chebyshev’s 
polynomials series in which coefficients are expressed through 
M l k (x) functions only (Ishov et al., 1989). Also it has been shown 
that all Lorentz functions M^(x) for k>2 can be limited to calcula- 
tions of the two first Lorentz functions M^(x) and M*j (x) only 
(Ishov et al., 1989). Moreover, the method of the universal func- 
tions has been developed for rapid estimating of the principal term 
of multil ine representation (4) (Shved et al . , 1 984, 1 988; Ishov et al . , 
1988). This method is based on introducing special absorption 
functions for each type of molecules (linear molecule and spherical, 
symmetric and asymmetric tops). For all these functions power 
(x *0) as well as asymptotic (X -»«>) series expansions managed to 
be obtained (Ishov, 1982). It was also shown that the universal 
functions for the Voigt line profile may be reduced to calculating 
corresponding functions for the Ixirentz and Doppler profiles 
(Ishov, 1985). 

I hus we can argue that high-precision and effective calculation 
of the principal term of multiline representation (4) is a solved 
problem. Although further study of mathematical properties of the 
absorption function Mq(x) for the Voigt line profile is, certainly, of 
applied importance. The thing is that different asymptotic represen- 
tations of M^(x) have irregular character with respect to the 
parameter a and do not cover the whole range of a and T parameters 
with sufficiently accurate calculations of Mq(x). 


MV(x) 

Mux) 
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Fig. 1 Isolines of relative error of the approximations used with 
calculation of the equivalent width of the band consisting of 

n=2, oo identical equidistant overlapping Lorentz lines: 

- the isolated spectral line approximation, 8j=| W N (n)- 
W(n)(l)| AV N (n); ... - the approximate version (6) of 
multiline representation (4), S 2 ==l W(n) ( 1 ) | /W N (n)+A(6) | / 
W N (n); - 8 3 “A(6)/W N (n). Numbers are corresponding 
to: 1-0.1%, 2-1%, 3-2%, 4-5%, 5-10%. Approximate repre- 
sentation (6) becomes exact in the points of intersection of 
and S 3 isolines. 8 is ratio of the distance between 
neighbouring line centers to the Lorentz half-width. 

5. NEW REPRESENTATION OF M^(t) 

Multiline representation (4) has also proved to be useful for 
deriving new representations of the absorption function Mq(t) for 
the Voigt line profile. Indeed, using the representation a v — 
(l-£)a D +£a L (i.e. assuming the each Voigt line consists of two 
lines and one with the Doppler profile and another with the Lorentz 
one) it may be shown that the following approximate relation is 
valid: 


Mo(i) - (1 - ^)qMg(t) + (x) + 

+ t[l-q(l-^ + ^)] > 

where q=H(0,V2a)/H(0,a), £=V2(l-q)/(V2-l) and H(x,a) is Voigt 
function. Relation (7) becomes exact as a=0 (Doppler profile case) 
and at a=°° (Lorentz a profile case) and has no peculiarities over the 
whole range of the a and z parameters. Relation (7) is worst satisfied 
at a= 1 , but even in this case it approximates the absorption function 
Mq(t) much better then the two first terms of Teylor’s series 
expansion. 

6. OTHER PROPERTIES OF (4) 

The multiline representation can be easily generalized for an 
optically inhomogeneous gas media: X-a(v^-v) — ~+J*dzSj 
(z)a(v^-v;z)n(z), where 1 is geometrical path length, n(z) is 
absorbing molecule concentration at the point z. With traditional 
application of the line-by-line method, the multiline representation 
eliminates exact accouting for the contribution of weak lines to total 
absorptance due to dividing all the lines into two sub-bands: strong 
and weak nes. The multiline representation holds true for combina- 
tions of molecular bands. 


7. CONCLUDING REMARKS 

Ihe present paper has revealed some basic properties of 
multiline representation (4) only. When calculating total absorp- 
tance in some spectral intervals covering spectral lines of several 
bands of various gases in different layers of the atmosphere, an 
alternative can arise within the framework of the proposed ap- 
proach. Particular attention should be given to approximate version 
(6) for multiline representation (4). The formulas (6) are simple and 
feasible for calculations. Therefore, it is desirable that the band be 
dividing into sub-bands so that each of them hit in the “ 0011 ^ 0 ^ of 
values 8 and T shown in Fig. 1, in w hich approximate relations (6) 
becomes extremely accurate. Besides, when writing approximate 
relations (5), the weakest line may be preferred to the strongest. In 
this case relations (6) will have the same form but be satisfied more 
exactly. 
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